Quantum vacuum radiation spectra from a semiconductor microcavity 
with a time-modulated vacuum Rabi frequency 
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We develop a general theory of the quantum vacuum radiation generated by an arbitrary time- 
modulation of the vacuum Rabi frequency for an intersubband transition of a doped quantum well 
system embedded in a semiconductor planar microcavity. Both non-radiative and radiative losses 
are included within an input-output quantum Langevin framework. The intensity and the main 
spectral signatures of the extra-cavity emission are characterized as a function of the amplitude and 
the frequency of the vacuum modulation. A significant amount of photon pairs, which can largely 
exceed the black-body radiation in the mid and far infrared, is shown to be produced with realistic 
parameters. For a large amplitude resonant modulation, a parametric oscillation regime can be also 
achieved. 
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The radiation generated by a time-modulation of the 
quantum vacuum is a very general and fascinating phe- 
nomenon, and has been predicted to occur in a variety 
of physical systems ranging from non-uniformly accel- 
erated boundaries (dynamical Casimir effect [3, 0]) to 
semiconductors with rapidly changing dielectric proper- 
ties [3[. These quantum vacuum phenomena have some 
analogies with the Unruh-Hawking radiation Q in the 
curved space-time around a black hole. Recent years 
have seen the appearance of a number of proposals to 
enhance the intensity of the quantum vacuum radiation, 
exploting, e.g., high-finesse Fabry-Perot resonators H, 
a time-modulation of the dielectric constant of a mono- 
lithic solid-state cavity [|| , or even the reflectivity change 
of semiconductor mirrors induced by an ultrafast photo- 
generation of carriers Q ■ Still, the very weak intensity of 
the emitted radiation has so far hindered its experimental 
observation. 

Planar semiconductor microcavities embedding a 
doped multiple quantum well structure have attracted 
a considerable interest in the last few years. As demon- 
strated by several spectroscopic experiments in the mid 
infrared range [1, 0, [HI U, 12, 3], the strong coupling 
between the cavity mode and the electronic transition 
between the first two quantum well subbands results in 
an elementary excitation spectrum consisting of inter- 
subband polaritons, i.e. linear superpositions of cavity 
photon and intersubband excitation states. The most 
interesting property of these systems from the point of 
view of the quantum vacuum radiation is given by the 
large vacuum Rabi frequency Q,r, being as high as a 
significant fraction of the intersubband transition fre- 
quency a; 12 11411 . In this unusual ultrastrong coupling 
regime 14J, [l5(, the antiresonant terms of the vacuum 
Rabi coupling play in fact a significant role and the quan- 



tum ground state of the system is a squeezed vacuum 
state containing a significant amount of correlated pairs 
of cavity photons. The photon pairs in the ground state 
are virtual and cannot escape the cavity if its parameters 
are time-independent (lH] . 

In order to release these bound photons into extra- 
cavity radiation, the quantum vacuum has to be mod- 
ulated in time. The recent experimental demonstration 
of a wide tunability of the cavity parameters (in particu- 
lar of the vacuum Rabi frequency ttji [Toj ) via an exter- 
nal gate bias and the possibility of ultrafast modulation 
[l2| makes the present system a very promising one in 
view of the observation of quantum vacuum radiation. 



A first theoretical study for an ideal isolated cavity [14 1 
has indeed suggested that a non-adiabatic switch-off of 
Qr results in a significant number of photon pairs being 
emitted from the cavity. 

A general theory which includes the effect of 
non-radiative and radiative losses for arbitrary time- 
modulations is however still missing, as well as a com- 
plete characterization of the intensity and spectral sig- 
natures of the emitted radiation for realistic cavity sys- 
tems. In this letter, we address these fundamental issues 
by developing a theory of the quantum vacuum radiation 
for an arbitrary modulation of the microcavity proper- 
ties. The properties of the extracavity emission for the 
most promising case of a periodic modulation of are 
calculated by means of the generalized input-output for- 
malism 15|: remarkably, the emitted quantum vacuum 



radiation turns out to be much stronger than the radi- 
ation by spurious effects such as black body emission. 
The instability regions in which the vacuum modulation 
produces a parametric oscillation of the cavity field are 
identified and shown to be within experimental reach. 

A theoretical description of the system can be obtained 
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FIG. 1: (color online). Top: a sketch of the considered semi- 
conductor planar microcavity system. Bottom: a scheme of 
the quantum model. 



by means of the formalism developed in [lj, [15J- The 
photon mode in the planar microcavity and the bright 
intersubband electronic excitation of the doped quantum 
well system (see Fig. [T]) are described as two bosonic 
modes. Given the translational symmetry of the sys- 
tem along the cavity plane, the in-plane wavevector k 
is a good quantum number. The creation operators for 
respectively a cavity photon and an electronic excita- 
tion of wavevector k are denoted by a k and S k . The 
in-plane dispersion relation of the cavity-photon is de- 
fined as w Cj k, while the frequency W12 of the intersub- 
band excitation is taken dispersionless. As explained 
in detail in Ref. the electric-dipole coupling be- 

tween one cavity photon and one bright intersubband 
excitation is quantified by the vacuum Rabi frequency 

1/2 
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mo "I3r ciei Ng f ^r fi 2 sin 2 6>J ' , where L e £ v 
is the effective length of the cavity mode, the dielec- 
tric constant of the cavity spacer, a e i the density of the 
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two-dimensional electron gas, i*q W 
ber of quantum wells, /12 the oscillator strength of the 
intersubband transition and 9 is the intracavity photon 
propagation angle such that sin 9 — c/c/(wi 2v / e c 7)- The 
second quantization Hamiltonian of the present cavity 
system reads 

lvivM k v k (1) 



H 



where the column vector of operators Uk is defined as 
«k = (Sk) ^k, «Lk> ^Lk) T ' V i s the diagonal metric 77 = 
diagfl, 1, — 1, — 1] and the Hopfield-Bogoliubov matrix 
Mk is defined as 
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For a quantum well, the diamagnetic coupling constant 
is approximately Z?k — ^Rk/^ 12 — !• '^' ne ultra-strong 



coupling regime is characterized by a value of f2ij jk 
comparable to w C)k and u>i2. In this regime, a cen- 
tral role is played by the anti-resonant light-matter cou- 
pling terms corresponding to the off-diagonal (1,3), (1,4), 
(2,3), (2,4) terms of M (and their transposed). In the 
following, a general time-dependence of the Rabi fre- 
quency is considered: Q R ,u(t) = ft RM + n™°*(t) and 
D k (t) =D k + D™° d (t). 

Non-radiative as well as radiative losses will be taken 
into account by means of the generalized input-output 



formalism developed in |15[: the system is in interac 



tion with two baths of harmonic oscillators, which are 
responsible for dissipation and fluctuations of both the 
cavity-photon and the electronic polarization fields. The 
radiative and non-radiative complex damping rates are 
denoted by_T Cik (u;) and ri2,k(^)- The real part (zero 
for lo < [15|) quantifies the frequency-dependent losses, 
while the imaginary part represents the Lamb-shift of 
the mode due to the coupling to the external bath. The 
resulting quantum Langevin equations are conveniently 
written in frequency space as the vector equation: 
00 

dcu' [M KiU 5{lo - J) + M£° d _ u ,] Skfy) = -iA 

O 

(2) 

where v^(u>) is the Fourier transform of the operator vec- 
tor Vk(t) and the quantum Langevin operator vector 



= (fc,kM,^kH,^_ k N,FL k (^) 

takes into account the quantum fluctuations due to the 
coupling to the baths. The matrix 

Mu,u = M k -lw 

- 1 diag[f Cik (o,), f 12 ,kM, f ;_ k (-o;), f ; a ,_ k (-a;)] 

summarizes the time-independent properties of the sys- 
tem, while -M k l ° d describes the time-modulation. For the 
case of a time-dependent vacuum Rabi frequency fi.R )k (t) , 
this has the form: 
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where ^^ k d w and D^ l ^f are the Fourier transforms of 
respectively f^° k d (i) and D™ od (t). 
The exact solution of |2]) is given by 



Vk(v) = ~i dui'g k (ui,ui')A(ui') , 



(3) 



where Gk(u>,uj') is the inverse of .M k (a>,u/) 



dw' V gZ°(cJ, uj')M!?{uj', u") = 5 rt S(u; - uj"). 
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Using the input-output scheme 15|, we obtain the spec- 
tral density of emitted photons from the cavity S£ ut (w) 
as a function of the incident one Sfp(u>) and the quantum 



Langevin forces T^: 
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(4) 



If one is interested in the quantum vacuum radiation due 
to the time-modulation of the cavity parameters, a vac- 
uum state has to be considered for the input state, so that 
(S™(u)) = 0, and the fluctuating quantum Langevin 
forces acting on the cavity-photon and the electronic po- 
larization modes € {c, 12}) are such that: 



k.k' 



W^M) = 47r«(w-cy)%,Re[T i|k (w)]<S: 




(5) 



This implies that only the last term of (3J gives a finite 
contribution to the emitted radiation. After some alge- 
bra, we get : 
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dw' {|a k 3 (w,-o;')| 2 K[f c ,k(^) 
|g k 4 (c,-c')| 2 ^[f 12 , k (c')]} ■ (6) 



The total number of emitted photons with in-plane wave- 
vector k is given by N° ut = f^S^iuj^doj. Note that 
in the absence of anti-resonant couplings in the Hamil- 
tonian (TTJ), C/ k 3 = C/ k 4 = giving a vanishing emitted 
radiation. 

This general theory can be applied to calculate the 
intensity of quantum vacuum radiation emitted by the 
cavity for an arbitrary modulation of the cavity param- 
eters and for arbitrary frequency-dependent losses. In 
the following, we shall focus ourselves on the case of a 
periodic modulation of the vacuum Rabi frequency, i.e. 



nrk(*) = AQ R , k cos(w mod t). 



(7) 



If the modulation frequency is tuned on resonance with 
the cavity modes, one expects [j| that the quantum vac- 
uum radiation can be strongly enhanced as compared 
to the case of a single sudden change of S!_R ; k discussed 
In the stationary state, the relevant quantity 
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to characterize the intensity of the emission is the total 
number of photons emitted per unit time dN£ ut jdt. 



Predictions for the rate dN£ ut /dt (in units of W12) 



of 



emitted photons as a function of the modulation fre- 
quency Lo m od are shown in the top panel of FigO For 



FIG. 2: Top panel: rate of emitted photons dN£ ut /dt (in units 
of W12) as a function of the normalized modulation frequency 
Parameters: (i J j C; - li + 2D k )/uJi2 = 1, T/lo^ = 0.025, 
fifl,k/wi2 = 0.2, Arijj ik /^i2 = 0.04. Note that, due to the 
scaling properties of the present model, the results do not de- 
pend on the specific value of W12. The letters A,B,C indicate 
three different resonantly enhanced processes. Bottom pan- 
els: the spectral density (arb. u.) for the processes A, B, 
C respectively. The resonant peaks occur at the LP (Lower 
Polariton) and/or UP (Upper Polariton) frequency. 



the sake of simplicity, a frequency-independent damping 
rate has been considered $t{T c ^(oj > 0)} = 5t{ri 2 ,k(^ > 
0)} = T; the imaginary part has been consistently de- 
termined via the Kramers-Kronig relatio ns II 51) . Values 
inspired from recent experiments @, [Tot Il2{ have been 
used for the cavity parameters. The structures in the 
integrated spectrum shown in the top panel of Fig[3]can 
be identified as resonance peaks when the modulation 
is phase-matched. Indeed, the vacuum modulation in- 
duces the creation of pairs of real cavity polaritons. This 
sort of nonlinear parametric process is enhanced when 
the phase-matching condition ru> mo d — w j,k + is 
fulfilled, r being a generic positive integer number, and 
j, f £ {LP,UP}. The dominant features A,B,C are the 
three lowest-order r = 1 peaks corresponding to the pro- 
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FIG. 3: (color online). Top panel: rate dN£ ut /dt (in units of 
W12) of emitted photons as a function of uJmod/u-12 for different 
values of the damping T/0J12 = 0.025 (solid), 0.05 (dashed), 
0.075 (dot-dashed). Other parameters as in Fig. [2] Bottom 
panel: normalized rate of emitted photons from a black-body 
emitter as a function of U12 for different temperatures. 
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FIG. 4: (color online). Top panel: rate dN£ ut /dt (in units of 
W12) of emitted photons as a function of w mo d/<^i2 for different 
values of the normalized modulation amplitude Af2_R i fc/cJi2 = 
0.01,0.04,0.07,0.1 (from bottom to top). Other parame- 
ters as in Fig. Bottom panel: instability boundaries for 
r/o)i2 = 0.025 (solid), 0.05 (dashed). Above the lines, the 
system is parametrically unstable. 



cesses where either two Lower Polaritons (LPs), or one 
LP and one Upper Polariton (UP), or two UP's are gen- 
erated. This interpretation is supported by the spec- 
tral densities plotted in the three lower panels of Fig[3] 
for modulation frequencies corresponding to respectively 
A,B,C peaks. In each case, the emission is strongly 
peaked at the frequencies of the final polariton states in- 
volved in the process; for the parameter chosen, we have 
indeed 1J, [l5| u>zp,k — ^12 — &rm = 0.8 ^12 for the 
lower polariton and u>up,k — ^12 + &R,~k = 1-2 W12 for 
the upper polariton. The shoulder and the smaller peak 
at LOmod/ ^12 < 1 can be attributed to r = 2 processes, 
while higher order processes require a weaker damping to 
be visible. 

More insight into the properties of the quantum vac- 
uum emission are given in Fig[3] In the top panel, the 
robustness of the emission has been verified for increasing 
values of the damping rate T, the resonant enhancement 
is quenched, but the main features remain unaffected 
even for rather large damping rates. In the bottom panel, 
comparison with the black body emission in the absence 
of any modulation is made: the total rate of emitted 
black body photons at given k is shown as a function 
of CJ12 (ranging from the terahertz to the mid infrared 
range) for k correspoding to an intracavity photon propa- 
gation angle of 60° and different temperatures. Note how 
the black body emission decreases almost exponentially 
with W12, while the quantum vacuum radiation, being 
a function of 0,R t u./^>i2 only, linearly increases with W12 
at fixed CIr^/wu- From this plot, one is quantitatively 
reassured that for reasonably low temperatures the quan- 
tum vacuum radiation can exceed the black-body emis- 
sion even by several orders of magnitude. The increase 



of the emitted intensity versus the modulation amplitude 
A^fl,k/o;i2 is shown in the top panel of Figf4] In particu- 
lar, note the strongly superlinear increase of the emission 
intensity around the A and C resonance peaks. In these 
regions, if the modulation amplitude is large enough, the 
system can even develop a parametric instability, the in- 
coherent quantum vacuum radiation being replaced by a 
coherent parametric oscillation [l6j]. Above the instabil- 
ity threshold, the results obtained from the solution of 
Eq. ^ in Fourier space are no longer valid, being the 
field amplitudes exponentially growing with time. Hence 
they are not shown here. The instability boundaries for 
parametric oscillation can be calculated applying the Flo- 
quet method 17j to the mean-field equations for the in- 
tracavity fields (ak) and (6k)- The result is shown in 
the bottom panel of Figf?]as a function of u> mo d/ V12 and 
Af2fl.fe/wi2: agreement with the position of the vertical 
asymptotes of the spectra in the top panel of of Figf4] is 
found. 

In conclusion, we have presented a complete theory 
with exact solutions for the quantum vacuum radiation 
from a semiconductor microcavity with a time-modulated 
vacuum Rabi frequency. In order to isolate it from spuri- 
ous effects such as black-body radiation, the main signa- 
tures of the quantum vacuum radiation as a function of 
the modulation parameters have been characterized. Our 
results show that semiconductor microcavities in the ul- 
trastrong coupling regime are a very promising system 
for the observation of quantum vacuum radiation phe- 
nomena. 
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